Fermion excitations of a tense brane black hole 
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By finding the spinor eigenvalues for a single deficit angle (d — 2)-sphere, we derive the radial 
potential for fermions on a d- dimensional black hole background that is embedded on a codimension 
two brane with conical singularity, where the deficit angle is related to the brane tension. From this 
we obtain the quasi-normal mode spectrum for bulk fermions on such a background. As a byproduct 
of our method, this also gives a rigorous proof for integer spin fields on the deficit 2-sphere. 
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I. INTRODUCTION 



o 

O, 

OO ■ Much has been said about black holes (BHs) in large extra-dimensional scenarios recently, e.g., see [l[ and the 

I references therein. Such theories lead to the intriguing possibility that the LHC could actually even produce BHs, 

. with various particle species being emitted by Hawking radiation or via classical BH excitations known as quasi- normal 

(~| modes (QNMs). Although it is not entirely clear as to whether or not detection of such processes is feasible, work still 

, remains to be done on various issues. For example, the effect of the brane tension for BHs with large extra dimensions 
has largely been ignored, because of the obvious difficulty of how to embed a BH onto a brane. 

O ' However, recently the work of Kaloper & Kiley inspired from codimension-two braneworld models, presented 

, the following metric for a black hole residing on a tcnsional 3-brane embedded in a six-dimensional spacetime: 

^ , ds^ = ~f{r)de + ^^+r'dnl , fir) = 1 - (^)' (1) 

where the radius of the horizon is given by 

^ : _ Mbh 



C — ■ . and Mbh is the mass of the black hole. The parameter 6 is a measure of the conical deviation from a perfect sphere 
■ and has the following angle element: 



dnj = del + sin^ 63 {del + sm^ 02 {d9j + b"^ sin^ eid<j)'^)) , 0<6<1 (3) 



^ , where for b — 1 this is the line element of the unit sphere S'^ and corresponds to zero brane tension. It may be worth 
mentioning that the location of the deficit angle is arbitrary and that it is possible to consider more than one deficit 
angle, where such cases may be of interest to the fermion generation puzzle re-expressed as a higher-dimensional 
problem, see Q and the references therein. 

For a non-vanishing brane tension the parameter 6 < 1 is a measure of the deficit angle about an axis parallel to 
the 3-brane in the angular direction (j>, such that the canonically normalized angle 4>' = (p/b runs over the interval 
[0, 27r/&]. It can be expressed in term of the brane tension A as: 
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where is the fundamental Planck constant of six-dimensional gravity. As can be seen the tension of the brane 
modifies the radius of the horizon, namely it increases with increasing tension (6 — *■ 0). 

The field equations for the scalar perturbations are discussed in [3, and similarly for graviton and electromag- 
netic perturbations in [gj. The emission rates were calculated in [3] for scalar, gauge boson and graviton fields, where 
numerical methods were used to solve the angular eigenvalue equations. Also, the scalar QNMs are found in [5[ by 
a perturbative expansion in powers of (1 ~ b) and hence are only valid for branes with small tensions. Recently the 
QNMs for scalar and gravitational perturbations have been found exactly, based on eigenflow arguments for integral 
1/b i. 

In this article we shall fill the gap by presenting results for the spin-half QNMs. Furthermore, the method we 
shall apply to spinors (by choosing azimuthal eigenvalue m = ±1/2, ±3/2, . . . ) also applies directly to the other 
perturbations (choosing eigenvalue m = 0,±l,±/2,...) and leads to a rigorous proof of the 2-sphere angular eigenvalue 
with no assumptions made on the form of b. 



II. SPINOR FIELDS 



The above works have only dealt with integer spin fields and in terms of phenomenology spin half fields are also 
important. For generality we shall begin our analysis by supposing a background metric which is d-dimensional and 
spherically symmetric, as given by: 



^fir)dt' + hir)dr'+r'dni^^, 



(5) 



where now dri^_2 denotes the metric for the (d — 2)-dimensional deficit sphere which in six dimensions has line element 
given by equation ([3]). 

Then under a conformal transformation: 



— _ o2 



where we shall take = 1/r, the metric becomes: 



dT 



h 



^dt^ + ^dr^ + dn^. 



■2i 



where 



tp — r 



(d-l)/2 



(6) 
(7) 
(8) 



(9) 



Since the t — r part and the {d — 2)-sphere part of the metric are completely separated, one can write the Dirac 
equation in the form: 



[(f Vt+rv.) ® 1] V' + 



^5 ^ f-a 



Sd-2 



V' = o, 



(10) 



where (7^)^ = 1. Note that from this point on we shall change our notation by omitting the bars. 

The problem now is to find the eigenvalue for the projected x\^\ which are the eigenspinors for the deficit (d — 2)- 
sphere, that is: 



(11) 



The separation follows exactly as in |i7|] and leads to the following radial Schrodinger-like equation in the tortoise 
coordinate r*: 



where G is the upper component of the two-component spinor 7] and the potential is given by: 



(12) 



(13) 



where we have set h{r) = l/f{r) for the Schwarzschild case, and the eigenvalue, k, shall be determined in the next 
section. 
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III. DEFICIT 2-SPHERE 



We shall first consider the case of a deficit 2-sphere, whereby we can generate results for the {d — 2)-sphere, using 
eigenflow arguments much like in 6]. First of all let's suppose the metric of the deficit two sphere be 

ds^ = dd^ + sin^ e d(t>^ , (14) 

where 6 is a positive real number and b — 1 represents a regular two sphere. The Dirac operator is then given by 



where the spin connection is given in terms of the zweibein and its inverse, 

where F" is the ChristofFel symbol. For the above metric, the only non- vanishing F" are, 



4"P 



-6^ sin 61 cos 61 ; FL cot 61. 



Choosing the zweibein to be 



and the Dirac matrices. 



= diag(l, 6sin6') ; ea'^ = diag(l, l/6sin6'). 



/=ai ; 7^ = a^ 



the spin connection are found to be 



r«=0 ; F, 



--bcos0 a^. 
2 



(15) 

(16) 

(17) 
(18) 
(19) 

(20) 



Here cr* are the Pauli matrices. 

Now the Dirac operator can be written down explicitly as. 



b sm 



2 / osmt' 



V'- 



Suppose we write the eigenvalue of this operator as zLin and express the fermion field ip in two component form: 



V'- 



Then we find the following set of equations: 



1 



6 sin 6* 



ip± = ±iK'4'±. 



Let us consider ip+, where ip- can be dealt with analogously. Consider the equation for 5^: 
that is, 

^(±) _ ^±im<p^ 



(21) 
(22) 

(23) 

(24) 
(25) 



Note that for spinors, one should get a sign change for a 27r rotation in ip. Therefore, the eigenvalues of m should be 
half-integers. 



1 3 5 



m 



2' 2' 2' 



(26) 
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At this point it may be worth mentioning that we can also obtain the scalar angular eigenvalue for the 2-sphere by 
assuming that m takes only integer values, i.e., m = 0, 1, 2, . . . . 

Returning to the eigenvalue equation for ■0+ we can make the following spinor separation of variables ansatz: 



r+nrn 



(27) 



Putting this into the eigenvalue equation, we have the following set of equations, 



do + l cote) 4±) T T^Mf^ = ^>^Bi^\ 



bsm9 



de + l cote] ± = 



6sin0 



These can be turned into second order equations. For Ait\ we have. 



do 



1 



cote 



b sin 9 



de 



1 



■ cot ( 



h sin 9 



One can get rid of the first derivative term by defining 

= (sin6')"^/^ 

Then the equation for u is 



i-4(f + i; 

IGsin^ ' 



l-4(f , ,2 

— 7^ ^ ^ 

16 cos^ 



2 -.-^^^^ 2 

This equation is just the one for the Jacobi polynomial with the solution 



u = 0. 



u = sm — 
2 



where P,"'^\coii9) is the Jacobi polynomial with 

m 1 

a = h - 

b 2 

and the eigenvalue obtained is 



2 



Pi"^'3)(cos0), 



TO 1 



TO 1 

to) = n + — + - 



(28) 

(29) 
(30) 
(31) 

(32) 

(33) 
(34) 



j(«:/3) 



where n = 0,1,2,... and to = 1/2,3/2,5/2,.... To ensure convergence n must be an integer and thus P„ is a 
polynomial. 

For the regular two sphere, with 6=1, we see that 



A = ?i + m + — =rt' + l. 



(35) 



where in the second step we have defined n' = 0, 1, 2, . . . with the constraint to = ±1/2, ±3/2, . . . , ±(ri' + 1/2), which 
is the expected result [8|. Thus, we can express the deficit eigenvalue as 



1 



K[n', m) = n+ \ni\ + - + - - 1 = n' + |to| 7 - 1 



1 



1 



(36) 



where we have generalized to include positive and negative to. 

It is now straightforward to generalize to {d — 2)-dimensions, because we know the result when b — 1 and hence can 
use eigenflow arguments similar to [^.^ Hence we find after dropping the prime: 



K(d, to) = n + 



d- 2 



ItoI ( - - 1 



(37) 



For 6 = 1 the standard result for the regular (d — 2)-sphere is obtained 0, 111- 



^ We have also confirmed this by checking the eigenvalue for the deficit 3-sphere. 
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TABLE I: Massless bulk Dirac fundamental (p — 0) QNM frequencies (Re{E) > 0) for a d = 6 tensional BH plotted for various 
{b,m,n) with jj. — 2. 



b 


n = 


n = 1 


n = 2 


n — 3 






m = 1/2 






1 


0.79441-0.39649 i 


1.30313-0.38936 i 


1.77861-0.39134 i 


2.24148-0.39212 i 


0.9 


0.79563-0.38139 i 


1.28425-0.37603 i 


1.74227-0.37791 i 


2.18878-0.3786 i 


0.7 


0.80602-0.34805 i 


1.24916-0.34612 i 


1.66789-0.3477 i 


2.07762-0.34822 i 


0.5 


0.83723-0.30914 i 


1.22496-0.30991 i 


1.59608-0.31101 i 


1.96117-0.31133 i 


0.3 


0.92644-0.26089 i 


1.24267-0.26211 i 


1.55178-0.26253 i 


1.8582-0.26263 i 


0.1 


1.35902-0.18211 i 


1.57048-0.18211 i 


1.7816-0.18211 i 


1.9925-0.1821 i 


0.01 


5.02563-0.08452 i 


5.12323-0.08452 i 


5.22083-0.08452 i 


5.31843-0.08452 i 


0.001 


22.71727-0.03923 i 


22.76257-0.03923 i 


22.80787-0.03923 i 


22.85317-0.03923 i 




m = 3/2 


1 




1.30313-0.38936 i 


1.77861-0.39134 i 


2.24148-0.39212 i 


0.9 




1.33616-0.37626 i 


1.79225-0.37803 i 


2.23806-0.37864 i 


0.7 




1.43038-0.34695 i 


1.84413-0.34801 i 


2.25202-0.3483 i 


0.5 




1.59608-0.31101 i 


1.96117-0.31133 i 


2.32389-0.3114 i 


0.3 




1.96004-0.26264 i 


2.26503-0.26265 i 


2.56951-0.26264 i 


0.1 




3.46647-0.18209 i 


3.67689-0.18209 i 


3.88729-0.18209 i 


0.01 




14.78521-0.08452 i 


14.8828-0.08452 i 


14.9804-0.08452 i 


0.001 




68.01597-0.03923 i 


68.06127-0.03923 i 


68.10657-0.03923 i 



IV. QNMS USING THE IYER AND WILL METHOD 

To evaluate the QNM frequencies we adopt the WKB approximation developed by Iyer and Will 9], also see 
references therein. Note that this analytic method has been used extensively in various BH cases [10|, where 
comparisons with other numerical results have been found to be accurate up to around 1% for both the real and the 
imaginary parts of the frequencies for low-lying modes with p < n (where p is the mode number and n is the spinor 
angular momentum quantum number). Furthermore, we have also included the p = n = modes in our results, 
shown in Table [J In figure [U we have plotted the fundamental QNM as a function of the tension, b. 

The formula for the complex quasi-normal mode frequencies E in the WKB approximation, carried to third order 
beyond the eikonal approximation, is given by [9']: 



= [Vo + (-2F„'V/2A] - z{p + i)(-2T/o'V/2(i + n), 



(38) 



where we denote Vq as the maximum of Vi and 



A = 



n = 





1 



288 \ Vr 



-^1 (77-fl88a2)--i^ 



Vr 



384 




(5)~ 



lOOa^) 



(19-f-28a^) - — 



(6)- 



288 \ Vn 



(39) 



(5 + 4.a'^) 



(40) 



Here 



a 



0,1,2, - •• , Re{E) > 
-1,-2,-3, ••• , Re{E) < 



and 



Vr 



(") 



(41) 
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FIG. 1: Plot of the fundamental QNM {p = 0, n — 0, m = 1/2) for varying tension, b. We see that as the tension increases the 
imaginary part vanishes. 



V. LARGE ANGULAR MOMENTUM & LARGE TENSION LIMIT 



A useful check of our numerical results is the exact analytic expression that can be obtained in the large overtone 
limit (also see @ for a similar discussion). If we now focus on the large angular momentum limit (k — > oo) we can 
easily extract an analytic expression for the QNMs to first order: 



Vo-t{p+^){-2Vo'y/^ + ..., 



(42) 



where Vq is the maximum of the potential Vi, see equation (|13|) . In the large n limit this potential now takes the 
form: 



where the maximum of the potential in such a limit is then found to be: 

'd-r 



Vn 



In this case wc find from the 1st order WKB approximation that 



E 



{d- l)2(d-3) th 



K — i{p+ -)Vd — 3 



(43) 



(44) 



(45) 



where this result agrees with the our previous result when b — 1 [7].^ 

Interestingly, choosing d = 6, and taking the limit b ~* leads to an expression independent of n (assuming m 
which is the case for spin half fields in six dimensions) 



E 



(5)i 



(46) 



where we have used the fact that rn — {^/bY^^ with fi = 2. Thus, as 5 ^ we see that the imaginary part becomes 
negligible (for fixed overtone p) and this agrees with our plot in Figure 1. Similar analysis in [6] shows that this limit 
is independent of the type of perturbation exciting the BH (to lowest order in an inverse power series in k). 



^ Note this is similar to the scalar field result [ill , but not identical. 
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VI. CONCLUDING REMARKS 

We have investigated the effect of brane tension on the low lying massless Dirac QNMs of a BH localized on a 
tense 3-brane in six dimensions. Conformal techniques were used to separate the time-radial and angular parts of the 
Dirac field analogously to the methods used in the zero tension case We have found that the only difference in 
the calculation between the tense and tensionless case appears in the value of the angular eigenvalues, k, which now 
depends on the the amount of tension, b. By calculating the unphysical deficit 2-sphere eigenvalue, k, we were able 
to determine the physically relevant 4-sphere eigenvalue by induction. We then computed the low lying QNMs for 
various tension, b, using the 3rd order WKB approximation see Table HI As a byproduct of our method, this 
also gives a rigorous proof of the result found in 6] for integer spin fields on the deficit 2-sphere with no assumptions 
made on the form of b, when wc choose integer values of m. 

As can be seen from the plot of the fundamental mode in figure [T] a general feature of these QNMs is that the 
imaginary part disappears with increasing tension (6 — > 0). Furthermore, as the tension is increased the real part 
grows without bound. These results agree with the large tension limit that was calculated analytically in section [V] 
Because the imaginary part gets smaller for larger tensions the amount of energy available for other bulk processes 
diminishes, which for example has been demonstrated in [3| for bulk BH emission rates for integer spin fields. Note, 
brane-localized QNMs/emissions are not affected by the brane tension, for a fixed horizon. 

Interestingly, the QNMs in the large p asymptotic limit can be calculated with tension using the method by 
Andersson and Howls [14], who have combined the WKB formalism with the monodromy method of Motl and 
Neitzke , also see [l6l|. The calculation in the tense brane case follows exactly that of the tenseless case 0] . Thus 
in the large p limit the QNM is purely imaginary with value, Ep = —i2i:THP- However, since the temperature is 
related to the Schwarzschild radius through the equation Th = {d — 3)/47rr^f and the radius is altered by the tension 
via equation ([2]), any tension will reduce the absolute value of the QNM in the large p asymptotic limit. 

In conclusion it is worth emphasising that tension is a fundamental property of branes in higher dimensional brane 
world models, yet until recent times, understanding the phenomenological consequences of tension on BHs has not 
been possible. Such knowledge is essential in order to determine realistic BH production rates and observational 
signatures at the LHC. In a further work we will investigate the Hawking emission for bulk Dirac fields from a BH 
embedding in a codimension two tense brane in six dimensions to compare these results with other fields on such a 
background and with fermions in the tensionless case p7j . 
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